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ABSTRACT 



Context. Most Wolf-Rayet stars (WR) of WC9 sub-type exhibit a dusty circumstellar envelope, but it is still a matter of debate how 
dust can form in their harsh environment. In a few cases, a pinwheel-like structure of the dusty envelope has been detected. Therefore, 
it has been suggested that dust formation in all dusty WR stars might be linked to colliding winds in a binary system. 
Aims. We probed the innermost region of the circumstellar dust shell of the deeply embedded WR star WR 118. 
Methods. We carried out spectro-interferometric observations using the AMBER instrument of ESO's Very Large Telescope 
Interferometer in low-spectral resolution mode (R = 35). The ^T-band observations were obtained with three 1.8 m telescopes span- 
ning projected baselines between 9.2 and 40.1 m. 

Results. At high spatial frequencies, the AMBER visibilities exhibit a prominent lobe, indicating that the envelope contains one or 
several zones with a large local intensity gradient. The strong closure phase signal clearly shows that the circumstellar envelope of 
WR 118 can only be described by an asymmetric intensity distribution. We show that a pinwheel nebula seen at low inclination is 
consistent with the AMBER data. Its size was determined to be 13.9 + 1.1 mas. 

Conclusions. WR118 possibly harbors a pinwheel nebula, which suggests a binary nature of the system. According to our best 
model, the period of the system would be « 60 days (for d = 3 kpc), making WR 118 the shortest-period pinwheel nebula known so 
far. 

Key words. Stars: individual: WR 118, Stars: Wolf-Rayet, Stars: winds, outflows, Stars: circumstellar matter, Techniques: interfero- 
metric, Techniques: spectroscopic 



1. Introduction 



Since the work of lAUen et al.l (|1972|) . infrared photometric stud- 
ies of Wolf-Rayet (WR) stars have shown that many late-type 
WC stars have dust characteristics in addition to their typi- 
cal free-free wind emission. For instance. [Williams et al.l ( 1987) 
showed that 85% of the WC9 stars and 50% of the WC8 stars 
are surrounded by heated (T ~ 1300 K) ci rcumstellar amor- 
phous carbon dust (see also lChiar et al.l l2002). WR dust makers 
are rare but remarkable in terms of their absolute dust-formation 
rate, which can be as high as M = 10" 6 M yr _1 . 

Among the strongest infrared WR sources, some show 
a distinct periodic variability. Variable dust producers ap- 
pear to be eccentric binary systems (WR+OB) where episodic 
dust formation coincides with periastron passage (see, e.g., 
Williams & van der Huchtl 1 1992b . Although our knowledge of 
dusty WR stars has considerably increased in the recent past, 
we are still far from a complete understanding of how dust can 
form in the hostile environment of these hot stars. 

Using aperture-masking interferometry with the Keck tele- 
scope, it has been discovered that two supposedly single WRs 
with large IR excess, WR104 and WR98a, are close binary 
systems whose dust distribution is tracing the orbital motion of 
the binary components and forming a pattern which has been 



* The observations presented in this paper were obtained with the 
AMBER instrument of ESO's Very Large Telescope Interferometer 
(VLTI) as part of the Guaranteed Time Programme 079.D-0359(A) (PI: 
T. Driebe). 



called "pin wheel" nebulae (|Tuthill et al.l 1 19991: IMonnier et al.l 
Il999ll2007|) . Recently, the number of WR stars and, in particu- 
lar, dusty WCs discovered has considerably increa sed, especially 
in some open clusters towards the galatic center (Ivan der Huchtl 
120061: ICrowtheret all 120061: ITuthill et al.l l2006h . Among them, 
two ne w pinwheel nebula e were discovered in the Quintuplet 
cluster (ITuthill et al.ll2006h . 

Binarity and wind clumping are currently the two best mech- 
anisms to form dust in the circumstellar environment of WR 
stars. Observationally, distinguishing between these two cases 
is a challenging task. Direct spectr oscopic evidence for a bi- 
nary companion is difficult to reveal ( Williams & van der Hucht 
2000), but infrared long -baseline interferometry with its spatial 
capability offers a unique chance to resolve a central binary star 
and to determine the spatial distribution of dust created in the 
circumstellar environment of dust-enshrouded WR stars. 

Recently, the dust shell of the persistent dust-maker WR 118 
(spectral type WC9), the third-brightest WR star in the K band 
(K = 3.65 m ), could be r esolved with speckle interferometry 
bv lYudin etakl (1200 lb and IMonnier et all (120071) . The radiative 
transfer modeling of lVudin et al.l (|200l|) showed that grains can 
grow to sizes up to -0.6 /mi. They also found that the inner 
dust- shell rim diameter of WR118 is r m = 17 mas, temper- 
ature of the carbon-rich dust at the inner dust- shell boundary 
is 1750 ± 100 K, and dust formation rate is of the order of 
10" 7 M o yr _1 . U sing the Keck t elesco pe, the aperture-masking 
observations of IMonnier et al.l (l2007|) confirmed the results of 
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lYudin etal.1 d2001b bv measuring a Gaussian FWHM of 23 mas 
at 2.3//m. 

Here, we present the first infrared long-baseline interfer- 
ometric observations of WR 118 obtaine d with AMBER, the 
near-infrared beam-combiner instrument (iPetrov et al .112007b of 
ESO's Very Large Telescope Interferometer (VLTI), giving ac- 
cess to spatial resolution (« 10 mas) five times better than before. 
This allowed us to investigate the geometry of the innermost cir- 
cumstellar dust shell of WR 118. 



2. Observations and data processing 

WR 118 was observed during three nights with AMBER in low 
spectral resolution mode (R = 35) in the /, H, and K bands. 
The telescope configuration used was E0-G0-H0 (linear array 
with baseline lengths of 16/32/48 m). Each measurement con- 
sists of 5000 frames with an exposure time of 100 ms. We re- 
duced the da ta with version 2.2 of the AMBER software pack- 
age amdlib dTatulli et al. l l2007lf l Following previous AMBER 
data reduction experience (e.g., I Wei gelt et al.l 12007k we kept 
20% of the frames with the highest fringe signal-to-noise ra- 
tio and discarded frames with an optical-path difference offset 
larger than 10 yum. Due to limitations related to the weather con- 
ditions and the extreme reddening of the object, the quality of 
the /- and if -band data turned out to be too low for a scien- 
tific analysis and were discarded. Thus, this paper discusses only 
the 7^-band (1.95 - 2.35/im) observations of WR 118. Also, we 
could only extract visibilities from the shortest baseline out of 
the 2007 dataset. An overview of the AMBER observations is 
given in Table [T] and the corresponding UV coverage is shown 
in the upper inset of Fig. [T] 



indicating that the circumstellar envelope deviates from a point- 
symmetric configuration. Therefore, from this qualitative anal- 
ysis of the AMBER data, we can clearly rule out any centro- 
symmetric shape of the object; i.e., the circumstellar shell can- 
not be described by a simple circu lar or even ell iptical envelope, 
as was done, for instance, in lYudin et al.l (1200 ll) . 



3. Modeling 

3.1. Simple geometrical models 

To interpret the AMBER data, we first tried to fit them with sim- 
ple two-component geometrical models: one resolved compo- 
nent (Gaussian disk, uniform disk, or ring) to account for the 
overall visibility decrease towards larger spatial frequencies, and 
one spatially off- set component (resolved or not) to account for 
asymmetry, mainly probed by the non-zero closure phase. We as- 
sumed that the object's intensity distribution is not wavelength- 
dependent in order to take advantage of the larger UV coverage. 
We performed a global optimization of the fit by using a simu- 
lated annealin g algorithm, comp lemented with standard gradient 
descent (see Millo ur et al.l l2009. for a first use of these algo- 
rithms). 



Table 1. Log of the AMBER observations of WR 118 (top) 
and the corresponding calibrators (bottom). B p gives the range 
of projected baselines, and PA denotes the position angle of 
the observ ation. The calibrators ' uniform-disk diameters were 
taken from Richichi et al.l {2005) and the spectral types are from 
Simbad (http : //simbad.u-strasbg . fr). 

The calibrated AMBER visibilities and closure phases of 
WR 118 are shown in Fig. Q] as a function of spatial frequency 
q. The figure reveals a visibility decrease up to q = 60, fol- 
lowed by a "bump" around q = 80. This visibility shape is char- 
acteristic for an object intensity distribution that contains sharp 
edges. Examples of such intensity distributions are, e.g., a uni- 
form disk or a ring-like intensity distribution. In addition, the 
AMBER data clearly show non-zero, non-180 closure phases, 
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Date 


Time 


Name 




PA 


seeing 




[UTC] 




[m] 


[°] 


["] 


2007-10-08 


02:18:53 


WR118 


9.2 


74.3 


0.61 


2008-03-28 


08:43:43 


WR118 


13.4-40.2 


61.7 


0.83 


2008-03-31 


06:51:05 


WR118 


9.6-28.8 


42.6 


0.73 


2008-03-31 


07:09:33 


WR118 


10.3-32.0 


47.6 


0.86 


calibrators 


Date 


Time 


HD 


UD 


spectral 


seeing 




[UTC] 


number 


[mas] 


type 


n 


2007-10-08 


01:07:58 


175583 


1.02 + 0.01 


K2III 


0.66 


2008-03-28 


09:33:52 


175583 


1.02 ±0.01 


K2III 


0.60 


2008-03-31 


07:42:46 


175583 


1.02 ±0.01 


K2III 


0.81 


2008-03-31 


08:20:04 


165524 


1.11 ±0.01 


K3III 


0.80 




Sp. Freq largest base (cycles/arcsec) 

Fig. 1. AMBER visibilities (top, points with error bars, UV 
coverage in the inset) and closure phas es (bottom) of W R 118, 
compared to the best-fitting model. The Yudi n et all (1200 ll) data 
(light blue) are shown in the top panel. The best-fitting geometri- 
cal model made of two rings and an extended background com- 
ponent (40 mas FWHM, not sketched) is shown as black lines 
(model image in the inset of the bottom panel). 
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Parameter 


Description 


Value 


SI 


I st ring size 


24.8 ±3.1 mas 


F2 


2 nd ring flux fraction 


34 ± 7% 


R. A. 


2 nd ring R. A.-offset 


0.79 ± 1.7 mas 


Dec. 


2 nd ring Dec-offset 


-6.0 ±2.3 mas 


S2 


2 nd ring size 


13.2 ± l.Omas 


F3 


Background flux fraction 


46 ± 11% 



Table 2. Parameters of our best-fitting geometrical model com- 
posed of two rings and an extended background component. 
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The reduced x 2 of the solution was typically 50 depending 
on the models we tried to fit to the AMBER data. Therefore, 
all these models are probably an inappropriate representation of 
the measured visibilities and phases. Thus, in the next modeling 
step, we tentatively added a new, fully resolved component that 
dilutes the flux of the whole system. By including this additional 
fully resolved component, the reduced x 2 f° r eacn model was 
considerably decreased (down to about 20). The results of our 
fitting procedure can be summarized as follows: 

It is not possible to find a reasonable fit to the data with only 
a two-component model. A significant improvement of the fit is 
achieved by adding a fully resolved background component to 
the two-component models, accounting for « 50% of the total 
flux. Thus, WR 1 18, as seen by AMBER, is composed of at least 
a fully resolved structure plus a component with a more complex 
shape, at smaller scales. 

All the models have a typical size of ^2 mas, which is ap- 
pro ximately in agreemen t with the results of lYudin et al.l (1200 ll) 
and lMonnier et al.l (12007b . 

The models containing an extended component plus a point 
source locate the point source inside the extension of the re- 
solved component. If this point source was a star and the re- 
solved component a dust shell, then the putative companion star 
would be located inside the dust shell. Then, the companion star 
would have probably carved out the dust shell. Therefore, the hy- 
pothesis of having a continuous dust shell plus a well- separated 
companion can be discarded. 

Our best-fitting model, which is able to reproduce both the 
visibilities and closure phases, consists of a set of two rings, with 
the smaller one located inside the larger one and with a contact 
zone in the South (see Tableland Fig.Q]). Best-fitting models 
using more than two rings also show interweaved rings of in- 
creasing sizes. 

Next, we provide a physical interpretation to describe the 
WR118 AMBER data, taking into account the knowledge ac- 
quired from the simple geometrical models. 

3.2. Pinwheel model 

iTuthill et al] dl999h and lMonnier etaD (ll999Ll2007h showed that 
the dusty WR stars WR98a and WR 104 are surrounded by 
spiral-like dusty nebulae. In particular, they found that the dust 
plume around these stars follow an Archimed ian spiral, pro- 
jected onto the plane of sky. lHarries et al.l (|2004|) presented a ra- 
diative transfer model of a bow- shock cavity between the winds 
of the two stellar components and successfully compared it with 
the observations. 

The goal of this section is to demonstrate that a pinwheel 
nebula is consistent with all the features seen in the AMBER 
data of WR 1 18. To explain these features, we developed a geo- 
metrical model of a pinwheel nebula, based on the same simple 
geometrical models presented in the previous section. It is gen- 
erated by adding many shifted elliptical rings of increasing sizes 
representing the bow- shock cavities at different orbital phases 
of the binary star along an Archimedian spiral, and then sub- 
sequently projecting them onto the plane of sky, as illustrated 
in the top panel of Fig. [2] Such a pinwheel model can be fully 
described by the following eight parameters: 

- the extension S , and the number of turns T of the spiral, 

- the size of the last ring s, corresponding to the opening angle 
6 - arctan^/S] of the wind- wind collision zone, 

- the flux fraction / of the last ring relative to the central one; 
the ring flux in the spiral linearly decreases from 1 to /, 
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Fig. 2. Top: Sketch of our pinwheel nebula model. A series of 
elliptical rings (three of them are highlighted) following a pro- 
jected Archimedian spiral (dashed line) set it up. The parameters 
of the model are labeled in the sketch. Bottom: Dependence of 
the pinwheel model on its parameters. For the first five parame- 
ters, the same model is shown in the first row. For the two last, a 
model with an inclination of 60° is shown. 



- the three projection angles: i, cd and y, to project the pin- 
wheel on the plane of sky, 

- F, the fractional flux of a fully resolved background. 

To illustrate the meaning of the individual input parameters, 
in the bottom panel of Fig. [2 we show how the morphology of 
the pinwheel changes if a sp ecific parameter is c hanged. 

Following the results of Tuthill et al. (2008), who showed 
that for WR 104, the pinwheel flux sharply drops after exactly 
one turn around the binary star, we decided to fix T = 1 and 
/ = (see Table [3]). Therefore, for our pinwheel model, we use 
the same number of free parameters (i.e. six) as in Sect. 13. II 

The best-fitting pinwheel model parameters are summarized 
in Table [3l and the comparison with the AMBER data is shown 
in Fig. [3] The qualitative agreement of the fit is very similar to 
the two-ring model described in the previous section, and the 
reduced^ 2 is also similar (x 2 = 16). Nevertheless, such a pin- 
wheel model is much easier to understand from a physical point- 
of-view, for instance, in terms of a dust-emission signature in the 
wind- wind curved shock around a binary star. 

Some parameters are well-constrained by our fit, such as the 
total extension of the pinwheel S = 13.9 + 1.1 mas, the extension 
of the extreme ring s = 22.6 + 4.7 mas, or the background flux 
contribution F = 49 + 5%. On the other hand, other parameters 
like the projection angles onto the plane of sky and especially the 
inclination angle are poorely constrained. This may come from 
the tight range of position angles of the AMBER data. 
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Fig. 3. Same as Fig. [T] but showing the best-fitting pinwheel 
model (black solid lines). An image of our best-fitting pinwheel 
model is shown in the bottom panel's inset. 



Parameter 


Description 


Value 


S 


Total size 


13.9 ± 1.1 mas 


T 


Turns 


1 (fixed) 


s 


Ring- size at extremity 


22.6 + 4.7 mas 


f 


Flux at extremity 


(fixed) 


i 


Inclination angle 


3 + 26° 


CO 


Orientation on sky 


9.2 + 9.6° 


7 


Rotation angle 


0.0 + 9.6° 


F 


Background flux fraction 


49 + 5% 



Table 3. Parameter of our best-fitting pinwheel model (see also 
Fig. El). 

We also found good agreement b etween our models and both 
the AMBER and lYudin et al.l (|200ll) data sets by using a size for 
the extended component of « 40 mas (see Figs. [T] and [3]), which 
is fully resolved by AMBER. 

3.3. Spherical wind clumping model 

A second physical hypothesis to form dust in WR sta rs is wind 
clumping in single star winds (iCherchneff et al.ll2000l) . We sim- 
ulated a spherical clumped envelope by randomly putting un- 
resolved "clumps" along a Gaussian profile with « 30 mas 
FWHM. If we assume a small number of clumps (< 10), the 
model visibilities hardly match the observed ones, while some 
clump configurations match the observed closure phases. For 
many clumps 1000, hypothesis favored by, e.g. jLepine et al.l 
2000), the model visibilities more easily match the AMBER 
ones, while, on average, the model closure phases show values 
closer to zero (or 180°) than the observations. 

Thus, from this study we qualitatively conclude that spher- 
ical wind-clumping models difficultly account for both the 
AMBER visibilities and closure phases. 

4. Discussion and conclusion 

Our study showed that our WR 118 data can be described by a 
dusty pinwheel nebula model. Other types of models describe 



the data equally well, but spherical wind-clumping models can- 
not easily explain both the observed visibilities and closure 
phases. Therefore, we suggest that WR118 is a binary star em- 
bedded in a pin wheel-like dust nebula. 

Assuming the pinwheel hypothesis, we can measure 
the opening angle of the bow- shock of WR118: 6 = 
arctan^/S] = 58 + 7°. This is c omparable to 40° in the 
case of WR104 dTuthill et al.l 120081) . We also infer a period 
for the pu tative binary system in WR118, assuming, like in 
WR98a dMonnier et al.l 1 19991) . tidal circularization of the or- 
bit. Given a typical wind terminal velocity of a WC9 star of 
1200 kms" 1 (ICrowthed [2007b. a distance of ^3 kpc for WR 1 18 
(Ivan der Huchtll200ll) . and an extension of 13.9 mas, the free- 
flying time of a dust plume from the center of the system to its 
farthest elongation is ^60 days, corresponding to one full rota- 
tion of the binary in our model. 

Such a period is remarkab ly shorter than the peri- 
ods of 241.5 days for WR104 (Ffuthill et al.l 120081). 565 days 
for WR98a dMonnier etZI 1 1999b and 24.8 vrs for WR112 
dMarchenko et al.ll2002l) . This would make WR 1 1 8 the shortest- 
period pinwheel system known today. 

The short period found for WR 118 may also be compared to 
other known WR binaries. For instance , y 1 Velorum is a W C8+0 
star with a similar period, of 78 days (Mill our et al.l l2007). but it 
does not produce dust. WR118 has a later spectral type than y 2 
Velorum, which alters the chemistry of the wind. Alternatively, 
the distance of WR 118 might be severely underestimated or its 
terminal velocity might be slower than for a typical WC9 star. 

With the study presented in this paper, we demonstrated the 
potential of infrared long-baseline interfero metry to resolve pin- 
whee l nebulae around dusty WR stars. In Ivan der Huchtl (2001, 
2006), we can find 24 permanent dust-producing WR stars, of 
which five (not yet including WR 118) have been confirmed to 
be colliding-wind binaries through the direct detection of pin- 
wheel nebulae. Eig ht have near-infrared sizes stron gly in favor of 
pinwheel nebulae dTuthill et al.l 120061: iMonnier et aUl2007b and 
are just waiting to be confirmed. In total, there are 18 dusty WR 
stars still to be resolved by long-baseline interferometry in order 
to provide observational evidence that all WR permanent dust- 
maker are colliding winds binaries. 
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